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Abstract. We introduce a generalized definition of parton distribution functions 
(PDFs) for a more consistent all-order treatment of power corrections. We present a 
new set of modified DGLAP evolution equations for nuclear PDFs, and show that 
the resummed asA^/'^/Q^-type of leading nuclear size enhanced power corrections 
significantly slow down the growth of gluon density at small-x. We discuss the relation 
between the calculated power corrections and the saturation phenomena. 



1. Introduction 

RHIC has produced good evidence tliat a new state of liot and dense matter of quarks 
and gluons, the quark-gluon plasma (QGP), was formed in ultra-relativistic heavy ion 
collisions py. To extract useful information on QGP properties, we need to understand 
the dynamics of parton multiple scattering and to derive precise nuclear PDFs (nPDFs) 
to calibrate the production of hard probes. In this talk, we use lepton-hadron deep 
inelastic scattering (DIS) as an example to demonstrate a need to modify the definition 
of PDFs for all-order treatments of parton multiple scattering and power corrections 
[2]. We present a new set of parton evolution equations consistent to QCD factorization 
beyond leading power and show that one-loop a^A^/^/Q^-type of power corrections 
significantly slow down the growth of gluon density at small-x. We argue that power 
corrections are very important for understanding the transition from a region where 
leading power pQCD has been very successful to the regime of parton saturation where 
new and novel strong interaction phenomena emerge [3]. 



2. QCD factorization and parton distribution functions 

Under the one-photon approximation, DIS cross section is determined by the hadronic 
tensor, W^'^{P, q), of hadron momentum P and virtual photon momentum q, represented 
by the imaginary part of partonic scattering diagrams in Fig. [H When = — is 
much greater than the hadron size 1 / fm, diagrams in Fig. [1] are dominated by the phase 
space near k"^ ~ (1/fm)^ <ti due to the perturbative pinch singularity at /c^ ~ [1]. 
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Figure 1. Feynman diagram expansion of DIS hadronic tensor W^, 



By applying the operator product expansion, the hadronic tensor can be systematically 
expanded in a power series of l/Q"^ and perturbatively factorized as [HI E] , 

W^'^ixB, Q') = E ^2;2n (g^, ^) ® /2+2n({x.}, /i^) (1) 

where Xb = Q^/2P ■ q, ® represents the convolution of parton momentum fractions 
Xi with i = 1 + 2n, /2's are PDFs, /i's with i > 4 are twist-z parton correlation 
functions (PCFs), and Hi'^ are corresponding photon-parton scattering amplitudes with 
all perturbative collinear divergences removed. In Eq. ([1]) and below, we suppress parton 
flavor dependence. The Tif^ can be perturbatively computed order-by-order in powers of 
as by applying Eq. ([T]) to a partonic state. At leading power, tree-level quark scattering 
amplitude provides parton model formalism. Radiative corrections at 0{as) can be 
derived, for example, in Fig. [21 by calculating the 0{as) quark scattering amplitude 
on the left and 0{as) quark distribution on the right. QCD factorization assures that 




Figure 2. Sample calculation of leading power partonic tensor at NLO. 



perturbative collinear divergence of the quark scattering amplitude on the left is exactly 
canceled by the same divergence of quark distribution on the right, and Ti-2(q) is infrared 
safe. The virtuality of the loop momentum k of the scattering amplitude on the left 
is limited by the kinematics at OIQ"^). The k in quark distribution on the right leads 
to a UV divergence that is removed by imposing a counterterm at a factorization scale 
/i^ ~ From Fig. O all low mass states of the quark scattering amplitude at 

< /i^ are canceled by that of the quark distribution on the right, and '^^2(5) keeps 
only the short- distance scattering dynamics at O^Q"^). The /i^ dependence of '^^2(9) 
an immediate consequence of the definition of parton-level quark distribution in Fig. [21 
and leads to DGLAP evolution equations of the PDFs. 

3. Power corrections and modified parton evolution 

To expand the reach of pQCD and explore the transition region to nonperturbative 
physics, systematic and reliable calculations of power corrections in Eq. ([ll) are crucial. 
When Xb Xc = l/2mR ~ 0.1 with hadron mass m and radius R, the virtual photon 
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in DIS could cover whole hadron in longitudinal direction to enhance power corrections 
[1]. Power corrections from tree-level diagrams were recently computed and found to be 
important at small xb The 0{as) correction could be computed in a way similar 
to the leading power. For example, in Fig. [3], is infrared safe with all coUinear 

divergence canceled between the left and the right. However, unlike the leading power 




k + . . . - \xb, ^2) + 



Figure 3. Sample calculation of l/Q^ power corrections at NLO. 

case in Fig. Ej the /c-integration of the diagram for the quark distribution on the right is 
UV finite, and consequently, ^^jj^] contains the physics of high mass states > Q"^) 
that were not in the partonic scattering amplitude on the left. 

Physically, we want all short-distance coefficient functions in QCD factorization 
formulas, like H.^^ in Eq. ([1]), to contain only dynamics of partonic scattering at Oi^Q"^). 
We can achieve this by defining the PDFs in Fig. [3] to contain only mass states with 

< /i^ ~ Q^. Consequently, from this additional /x^ dependence of PDFs, DGLAP 
parton evolution equations receive power corrections [2], [8] 



_d_ 



/2(x,/i^ 




X 
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/2+2n({l/i},/i^) 
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where /j with i > A obey similar equations to form a closed set of evolution equations 
[2]. The P's can be perturbatively computed from /i^ dependence of PDFs and PCFs. 

If the parton momentum a; 3> with saturation scale Ql{xs) = /i^ [3], the one- 
loop contribution to quark distribution in Fig. [3] is dominated by diagrams with two un- 
pinched poles that fix two of the three parton momentum integrals. Under this leading 
pole approximation [9], only diagrams with maximum un-pinched poles contribute and 
all Hi but one in Eq. ([2]) are fixed by these poles. We computed all one-loop leading 
OsA^/^/Q^-type of power corrections to parton evolution of nPDFs [2] 

d ( q{x,iJ,^) \ ^ dy 
^ Ix y 
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2tt 
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where Ta is a translation operator that shifts momentum fraction y — » y{l + A) with 
A = gjA' = ^2(^1/3 _ iyg2 ^^^^ ^2 ^lefined in [7], and Pgg, Pqg, Pgg are normal leading 
order splitting functions, while Pgg is the normal g ^ g splitting function with the term 
proportional to Uf excluded. In Fig. 111(a) and (b), we show gluon evolution without 
(red) and with (blue) power corrections, and the ratio of these two different evolutions, 
respectively. Using three sets of available PDFs, we found that the resummed power 
corrections can give as large as 20% reduction in slope of gluon evolution and can 
significantly slow down the growth of gluon density at small- a; and low Q^. 
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Figure 4. Modification of gluon evolution for nucleus oi A= 197 



4. Summary and outlook 

We modified the definition of PDFs so that all short-distance coefficient functions of 
power corrections in QCD factorization formulas contain only dynamics of partonic 
scattering at OiQ"^). Consequently, we obtained a new set of parton evolution equations 
including calculable power corrections in Eq. ([2]). Under the leading pole approximation, 
we computed all one-loop contributions to the evolution equations of nPDFs and found 
that leading a^A^/^/Q^-type power corrections significantly slow down the growth of 
gluon density at small-x and drive the gluon density to saturation [2J. 

Our factorization treatment of power corrections can be applied to any factorizable 
observables. Power corrections closely connect to coherent multi-parton dynamics and 
are important for studying the transition from a dilute to a saturated partonic system. 
When partons are close to be saturated, or all Xi are close to the Xs, the leading pole 
approximation is no longer reliable. Many more diagrams will contribute. We need 
different approaches for studying power corrections and parton correlations [H [8] . 

We acknowledge support from the U.S. Department of Energy Grant No. DE- 
FG02-87ER40371 and Contract No. DE-AC02-98CH10886. 
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